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Investigations were carried out to analyze various properties of materials em- 
ploying computer simulation techniques based on semiempirical potential functions 
>vhich comprises two-body and three-body interactions, In the following three ‘self 
contained' chapters we outlined the basic procedures employed and the advances 
made during this period . 

The first chapter involves an analysis of the relationship between the thermo- 
dynamic stability and the parameters of the semiempirical potential energy function 
employed in the simulation calculations. It contains a detailed investigation of the 
structural transformations from Diamond Cubic => /34in and from Graphite 
Diamond Cubic which occur at high pressures. Also, in this chapter the surface 
reconstruction of the (lOO), (110) and (111) planes for the diamond cubic crystal 
was analyzed and a marked contraction of the first surface layer spacings was found. 
Furthermore, it was found that the reduced aurface tension, o”, decreases with in- 
crease of 2', the three-body intensity parameter, for all three planes and falls in 
the expected experimental range. 

The second chapter involves a parametric analysis of the elastic constants and 
the mechanical stability of cubic structures against all types of homogeneous elastic 
deformations. The domain of absolute stability, involving simultaneous thermo- 
dynamic and mechanical stability for fee, bcc, diamond cubic and simple cubic 
structures Were determined using a potential energy function based on tv/o- and 
three-body interactions. 

In the last chapter, specific parameters involved in two-body plus three-body 
semiempirical potential energy function for carbon, silicon and silicon carbide sys- 
tems were calculated. Excellent reproduction for various energetics and structural 
quantities were obtained. 

These three chapters were prepared as separate manuscripts and submitted for 
publication. 
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Chapter 1 


A STUDY OF CRYSTAL STABILITY AND SURFACE ENERGY 
FOR DIAMOND CUBIC STRUCTUTRES 
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IHTRODOCTION 


From intermolecular iiorcc studies, it lo now kn^wn that the overall 

non-additlve contribution to the lattice energy Is positive so that 

analysis based on only pairwise additivity suggests a shallower inter** 

molecular potential than the true value^^\ Two body contributions alone 

are also know to be categorically unable to even qualitatively describe 

some configurations of molecular clusters In the. gas phase^^^ or the 

general relaxation and reconstruction of fee crystal surfaces^^\ In 

addition, the many-body contribution has been shown to play a key role in 

the stability of certain crystal structures^^\ 

( 2 - 4 ) 

In these recent analyses' a relatively simple potential energy 

function (PEF), comprising only a two-body tile-type potential plus a 
three~body Axilrod-Teller-type potential, was found to be extremely 
effective. In the present paper, this same parametric PEF is applied to 
describe the bulk stability and surface energy for the diamond cubic 
structure. To test the stability condition, the FCC, BCC, diamond cubic, 
graphite and p-tln structures were considered. For the surface energy 
study, only the (IH), (100) and (110) planes of the diamond cubic 
structure were considered. 


3 


ANALYTICAL PROCEDOSE 
(4) 


The PEP, 1 ^, le given as 


$ - 4<2 *•* < 1'3 

where 


and 
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In eq. (lb), r^j »« | “ r^ ) , denotes the equilibrium atomic 
separation distances, e Is the two-body energy at r. . ■ r while the 

X j O 

exponents m and n account for the fepulslve and attractive two-body 
contributions respectively. In eq. (Ic), 0^, 0^, and r^^.. , r^j^, r^j^ 
represent the angles and sides formed by three particles 1, J and k 
while Z is the three body Intensity parameter. In this general form, 

Is a function of both the material parameters e, r^, m, n, Z and the 
atomic configuration of the system. 

Choosing m « 12 and n * 6, using unitless quantities, and lattice 
sums, the total reduced Interaction energy for a crystalline structure, can 
be written from eq. (1) as 
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with d denoting the nearest neighbor distance In the crystal. For each 
chosen crystal structure, the lattice suns have unique values which are 
listed In Table 1. 

The convergence of the two body lattice suns, A^, Is rapid for large 
m. Computed lattice sums Ag are In good agreement with the 

results of Bell and Zucker for the BCC, FCC and HOP structures^^\ The 
sun Tj^ is very slowly converging. 

The stability condition for a crystal at T *• 0°K can be obtained by 

A ^ 

considering 6? /5v "0 or 3® /3d ■ 0 since the atomic volume v Is 
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related to d by v " gd where g Is a geometrical constant. Values of 
g are tabulated In Table 1 for various crystals. Thus, the stability 
condition becomes 


Aj2<r*)^^ - Ag(r*)^ | 2* \(r*)® - 0 . (3) 

The surface energy, a, per surface atom was calculated from 
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a 


(4) 
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where is the total Interaction energy for an atom located Ip the i'th 

layer from the surface, Eq Is the corresponding value for a layer far 

inside the crystal away from any Influence due to the surface, M is the 

total number of surface layers considered in the calculation and N is 

a 

the total number of atoms at the expected surface. At a finite 
temperature, <s is expressed as the sum of two parts, a and 6 ,, 

where represents the surface energy for the atoms in unrelaxed 

positions while the change due to relaxation given by 


‘rel ■ n; - -“o’ • 

Throughout this investigation, relaxations for M' ■ 7, and T - 298“ K 
with the remainder of the atoms fixed were performed using a Monte Carlo 
technique based on the Metropolis approximation^ 


RESULTS AND DISCUSSION 

A. Structural Stability 

Using eqs* (2) and (3), minimum values of 9 , for various crystalline 

* 

structures, were calculated as a function of Z and are shown in Figure 
1. We note that the more densely packed structures are favored at small 

iV ^ 

Z and the more loosely packed structuires at large Z . The diamond cubic 

ift 

structure Is the most stable in the range 0.55 < Z < 0.80. Including the 
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(4) 

previous results of Hsllcloglu that the HCP structure has the lowest 
energy without the three body Interaction (l*e« > for Z 0), structural 
stability appears to change In the following order with Increasing Z* 

HCP FCC •»“ DIA + GRAPH » <6) 

ik 

we were unable to find any region of Z where the BCC structure Is the 

most stable structure for m <* 12 and, n ■* 6; however, the BCC structure 

becomes stable for other values of m and n« 

Figure 2 shows our calculated results for as a function of atomic 

volume for Z •• 0*3 (FCC most stable), Z "0.7 (OIA most stable and 

Z* " 11 (GRAPH most stable). We note that both the minimum energy and the 

A ^ 

equilibrium volume, v , (volume at Increase for all structures 

it * 

as Z Increases. Table 2 lists these numerical values of and 

* 

v for the different structures. Considerable variation of these values 
* 

with Z exists for the FCC structure but not for the graphite structure. 
From Fig. 2b we note that, for the stable diamond structure, If v is 
decreased by 14% or more (by Increase of pressure) a phase transition Is 
expected to occur to the (J-tin structure. Likewise, from Fig. 2c, one may 
anticipate a graphite •*’ diamond structural transition If v Is decreased 
by 19% or more. Thus, these parameterlc studies lead to the conclusion 
that, within the stability region of the diamond cubic or graphite 
structure, the diamond cubic + p-tln and the graphite diamond cubic 
structural transformations may be expected upon application of sufficient 
external pressure. 


/ 


% 


Most group XV materials (C, SI, Ge, Sn) exhibit the diamond cubic 

structure and their phase diagrams show that SI, Ge and Sn tend to form 

the p-tln structure under high pressure* Furthermore, the stable 

structure for C at NTP Is the graphite structure which tends to 

transform to the diamond cubic structure under sufficiently high pressure* 

By evaluation of the clastic constants, the diamond cubic form has 

* 

been shown to be absolutely stable for Z > 0.55 but unstable to shear 
forces below this value. Thus, mechanical stability of this crystal 
structure requires the presence of a many-body force* This study will be 
published elsewhere^^\ 

B* Surface Energy 

For the relaxation procedure, the surface and bulk energies, o and 

A 

Eq were calculated for different Z values ranging from 0*6 to 0*9 where 

the diamond cubic structure Is stable* The (111), (100) and (110) surface 

planes were chosen for consideration and the top-most two surface layers 

were relaxed to minimize o* The calculation system contained 

approximately 260 atoms and considered about 10 atomic layers. In Table 3, 

values of and a are given for these three orientations In 

reduced units at various Z*. We note that the (111) has the lowest 

surface energy while the (100) has the highest and that cr decreases for 

A 

all three orientations as Z Increases. 

The top and side views of the (100), (110) and (111) planes with 
Z* ■ 0.7 after relaxation are shown in Fig. 3. For all three surfaces, the 
top layer exhibits a large contraction while changes In the second layer 
relative to the bulk spacing are very small. The surface geometry after 
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relaxation waa found to be relatively insensitive to Z value in the 
diamond cubic stability range. No change with relaxation waa found for the 
top view of the (Ul) and (UO); however, the (100) exhibited dimer 
formation. Tliis reflects the relaxation energy, which Table 3 

shows us Is 5 - 10 times larger for the (100) than either of the other 

surfaces. The reconstruction patterns for the (100), (110) and (111) 
surfaces are thus C(2 x 2), (1x1) and (1 x 1) respectively. 

In Pig. 4, the reduced surface energy a* - Ict/EqI and the bulk 
energy. Eg, are plotted as a function of Z*. For all three orientations, 
a was found to decrease monotonically with Increase of Z- . The bulk 

A 

energy, Eg, exhibits the opposite trend with increasing Z . Experiment 
tally, approximate values of reduced surface energy for Group IV elements 
(C, Si, Ge and Sn) are in the range O.U - 0.25 and these are within 
the calculated range of Fig. 4. 

To estimate the "relative" experimental surface energy, we take the 

raw data of Table 4 for lattice constant, Sg, surface energy, Yt 

enthalpy of forrsatlon at 298®K, Since the surface area per 

2 2 2 

surface atom a is given by 0.217 Sg, 0.354 Sg and 0.5 Sg for the 
(111), (110) and (100) planes respectively, we can define the reduced 
surface energy by 9 “ ya/AH^^g. Values of 9 for the major orientations 
and the average value 9 are given in Table 5. 

COIiCUISIONS 

Atomistic calculations using a two-body-only PEF are unable to predict 
the relative stability for FCC, BCG, diamond cubic and graphite structures; 
however, addition of a many-body term in the form of a three-body 
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contribution does allow such a prediction, the present calculations 
« 

predict that the structural stability order with increasing three-body 

strength, S*, Is HCP FCC p-tln Diamond Graphite. In addition, 

structural transformations from diamond cubic ^ p-tin and graphite 

diamond cubic are predicted to occur by Increases of pressure. 

Multilayer surface relaxation of the principle orfcintatlons In the 

diamond cubic structure shows marked contraction of the first Interlayer 

spacing for the (100, (110) and (111) with dimer formation on the (100) and 

A 

calculated surface energies, a , In the expected range for the experimental 
values. 

This study emphasises the need for Inclusion of multibody forces In 
the PEP when dealing with both wtructural stability and surface energies of 
non-close packed structures. The PEP employed In this study was a 
particularly simple and useful one because It allowed ease and speed of 
computation In lengthy Iterative procedures. 
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Table 1 

Lattice Sums for the Mle and Axilrod-Teller Potentials 
Plus the Geometrical Constant, g, for Various Structures 



PCC 

BCC 

DIA 

GRAPH 

MB 

^6 

14.4481 

12.2495 

5.1153 

3.3895 

8.2864 

^12 

12.1319 

9.1141 

4.0389 

3.0092 

5.4654 

B 

19.1697 

14,7719 

1.6647 

0.1010 

7.0706 

If 

8 

0.70711 

0.7698 

1.5396 

3.067 

0.9797 


6912 

8192 

5832 

5040 

7200 


t indicates the number of atoms which were considered to 
compute Che lattice sums. 

fig ■ v/d^. 
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Table 2 

• 

The Dlnlmum energies, and the equilibrium volume, 

A t 

V , for five different etructures as a function of Z « 


* 

z 

★ , it 

PC 

BCC 

Diamond 

Graphite 

p-tln 

0.3 

®mln 

ff 

-3.8212 

-3.7494 

-2.6086 

-1.8731 

-3.4417 


V 

0.6925 

0.9051 

1.4854 

2.9105 

1.0057 

0.7 

9* 

'^min 

it 

-1.5049 

-1.S871 

-1.9763 

-1.8265 

-1.7145 


V 

1.3065 

1.3116 

1.6562 

2.9382 

1.3475 

1.1 

^mln 

ft 

-0.0073 

-0.8165 

-1.5188 

-1.7812 

-0.9640 


V 

1.7835 

1.7806 

1.8425 

2.9662 

1.7426 

t All 

values 

am reduced according to 

eq. (2). 
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Table 3 

CfAculated surface energies of Che low Index planes 
for Che diamond cubic scruccure aC several dlfferenC Z* values 



(100) 

(110) 

(111) 

<^un 

0.8266 

0.4374 

0.2084 

^rel 

'•0.3624 

-0.0298 

-0.014 

a 

0.4342 

0.4076 

0.1944 


0.7426 

0.3630 

0.1696 

rel 

-0.3714 

“0.0492 

-0.0292 

a 

0.3712 

0.3138 

0.1404 


0.6672 

0.2978 

0.1362 

^rel 

-0.3872 

-0.075 

-0.0406 

CT 

0.2800 

0.2228 

0.0956 


0,.6010 

0.2414 

0.107 

*rel 

-0.4152 

-0.1114 

-0.1006 

a 

0.1854 

0.130 

0.0065 


To cal number of 
atoms In Che system 

216 

360 

288 

Total number of 




atoms on the surface 

18 

30 

48 



0 

* 

Table 4 

>Jlrdup IV experimental data for lattice constant » a^, 
surface energy, y» and bulk energy, AH^gg 


Element 

ao<A) 

Y(ergs/cm ) 

4H298(eV/atora) 

n 

diamond 

3 . 5668 ^^^ 

3760 . 0 ^®^ 

7.412(10.11) 

SI 

5 . 4307 ^^^ 

1389 . 0 ^®^ 

4 g73<10,U) 

Ge 

5 , 6574 ^^^ 

975 . 5 ^®^ 

3 . 85 ^^^^ 

Sn 

6 . 4912 ^^^ 

822 . 0 ^®^ 

3.14<^2) 


15 




Table f 


Reduced surface tension, 6, for Group IV elements 


Species 

*(111)(cm^) 

*(110)(cm^ 

“(100) (cm^) 

®(111) 

®(110) 

‘'^*“(111) (eV) 

''^*®(110)(eV) 

‘'^'“(lOO) (eV) 

^diamond 

2.7544 X 10“’® 

4.4979 X 10“’^ 

6.3610 X 10“’® 

0.0873 

0.1426 

0.6472 

1.057 

1.4948 

SI 

6.3853 X 10“’^ 

1.0427 X 10"’* 

1.4746 X 10*’ * 

Oe 11862 

0.1937 

0.55432 

0.90519 

1.28014 

Ge 

6,9295 X 10“’® 

1.13157 X 10“’* i 

1.600 X 10"’* 

0e1097^^ 

0.17919 

■1 

0.6899 

0.9755 

Sn 

9.1226 X lo“’® 

1.4897 X 10“’* 

2.1068 X 10"’* 

0.1493 

0.2W 





®{100) 


1 

0.20167 

0.1439 

0.2739 

0.1954 

0.2534 

0.18077 


.3W 0.2459 


0.46B67 


0.76533 


1.08237 
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FIGURE C&FTIOHS 

it 

1. The minimum reduced energy, , for FCC, BCG, diamond cubic (DIA) and 
graphite (GRAH) structure as a function of the reduced three-body 
parameter, Z*, (m * 12, n ■ 6). 

it 

2. The reduced energy^ 4 » as a function of the atomic volume i V, for 
FCC, BCC, DIA, GRAH and p-tln structures (m - 12, n ■ 6); 

(a) for Z* - 0.3, within the stability region of FCC, 

(b) for t > 0.7, within the stability region for DIA, and 

(c) for Z* « 1.1, within the stability region for GRAH. 

3. Surface geometries (top and side views) of (a) (100), (b) (110), and 
(c) (111) surfaces for the diamond cubic structure after multilayer 

it 

relaxation with Z *0.7 and Monte Carlo temperature of 298°K<. 

The numbers below the figures Indicate the percent changes In 
Interlayer spacing. 

it 

4. Plot of the reduced surface energy, c , for three different surfaces 

and the bulk energy Eq on a per atom basis In the DIA system as a 
* 

function of Z at T ■ 298 K. 
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FIGURE 2a 
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Chapter II 


ABSOLUTE CRYSTAL STABILITY AND ELASTIC CONSTANTS 

FOR CUBIC STRUCTUTRES 


ItmODOCTION 


t 

In a pravlous paper, the relative atablllty of varioua open crystal 
structures was shown to require a multi-body potential interaction in 
addition to the usual two-body Interaction. There, the toultl-body 
term was chosen to be a simple three-body potential and a parametric 
analysis was used. In the present work, this analysis Is extended to 
examine the mechanical stability of various cubic structures against all 
types of homogeneous clastic deformations. As a necessary precursor, a 
parametric analysis of the elastic constants is obtained. 

Using unitless quantities, a two-body Mle potential, a three-body 
Axllrod-Tcller potential and lattice sums of the energy, the total reduced 

it 

interaction energy, ® , for on atom in a particular crystalline structure 
was shown to be given by^^^ 


* 

9 


+ Z* 


with 


9 




* 

Z 



and 


A 


N . „ K 

I (t^)" , b - J < 
j>l ’^ij J>1 


d ^n 






, N N 

4 1 1 

^ 1 i 
i>j>i 


(1+3 cos 6^ cos 


0. cos 0.) d 

^ i 






(la) 


(lb) 


(lc> 
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Horot Dr Hi c and are the usual two-'hody parameters while Z is the 
three**body Intenalty parameter and d la the nearest neighbor dlstonee in 
the crystol* For each chosen crystal structure, the lattice sums, A, B 
end in eqs. (1) hove unique values and their computed values are 
listed in Table 1 for the various cubic structures of interest to us in 
this paper* 

The equilibrium nearest neighbor distance, d^, at T •> 0**K may be 
found from eqs. (1) by requiring that the total energy be a minimum with 
roppeut to the nearest neighbor distance, d. As found earlisr^^^ the 

I 

equilibrium equation is then written as 

A(r*)“ - B(r*) ^ + 18 z* T. (r*)® - 0 . (2) 

Using eqs* (1) and (2) plus elasticity theory, the elastic constants, 
mechanical stability conditions and the absolute stability conditions for 
various crystal structures using this parametric potential may be obtained. 

ELASTIC CONSTANTS 

In any study of mechanical crystal stability the elastic constants are 
( 2 ) 

an essential quantity. ' They may be calculated by applying a 
homogeneous deformation to the crystal and comparing the stored; energy 
change with that predicted by linear continuum elasticity theory. In this 
procedure, the position of the I'th atom is denoted by in the final 

( 3 ) 

configuration. According to Wallace' ' ^ the change in distance between 
any two atoms in a crystal under homogeneous strain can be represented 
using Fig. 1 by 
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(3) 


I (xj)' - I uy + 2 q xj x' . 


where the indices a, ft represent Cartesian coordinates and take on the 
values X, y, 2 or 1, 2, 3» and the tij^p are the Lagranglan strain 
parameters. 

The elastic constants at 0°K are defined as 


r - i f ^ ^ ) 

<=afr6 V - 


where V denotes the volume of the crystal. These elastic constants are 
calculated at the initial configuration, with the strain T)^p being 
measured from r^ and the derivatives all being evaluated at r^, 
l.e., at Tijjp " 0* Using eqs. (1) and (4), the reduced elastic constants 

A 

may be written, using - (V/Ne) as 


* mn ym , , v . mn ^n , . „ 


* y *v9 

" ’•aPY6<'^ > • 


The lattice sums for the elastic constants, ^^pyp* ®opy6 "^oPy6’ 

are developed in Appendix X, have been computed and are listed In Table 2 
for various cubic structures. 

ik 

There are 81 elements in the matrix and, considering the 

(4) 

strain and crystalline symmetries plus using the Voigt notation , the 
matrix of clastic constants for the cubic system becomes 


27 


I 



In Eq. (6), there, are three Independent elastic constants; l*e, 

lit llf 

C^2 ^44* well known that the Cauchy relations, C^2 ” ^44 

for the cubic system, hold if each pair of atoms interact through a purely 
central potential^ . This relationship can be seen in Table 2 for the 
two-body part of the lattice sums for the elastic constants, i.e,; 

Aj 2 " and Bj2 “ ®44* course, it does not hold for the three-body 

portion which is not a central potential. 

Using eqs. (2) and (5) with m >■ 12 and n ■■ 6, the elastic constants 
C* for the FCC, BCC, diamond cubic and simple cubic crystal structures 

it 

have been calculated as a function of Z and are shown in Figs. 2. 

Except for the Cj^2 ^44 simple cubic structure, all the 

it 

elastic constants, C^2> ^44^ , decrease as Z Increases. 
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ABSOLDTB CKTSTAL ST^ILITY 




For a crystal to be stable against all elastic homogeneous deforma- 
tions, the elastic strain energy, must be positive. 

This condition requires that the matrix of elastic constants, (C ] , be 

ic 

positive definite; l.e. , all eigenvalues of the matrix 


(4) 

positive. Thus, the condition for mechanical stability In the cubic 
system Is 


C^i > 0; C^j>|Cj 2 l; C ^4 > 0; + 20^3 > 0 . (7) 

Calculated values of - 1^121 each structure 

have been plotted In Fig. 2 to assess the mechanical stability condition 
via eq. (7). 

A 

For the FCC structure, a value of Z < O.S satisfies the stability 

condition giving an extended domain of absolute stability. We note that 

this potential function gives no region of mechanical stability for the BCC 

structure although In another study' * It was found that altering m and 

n led to a stability domain for the BCC structure. We note also that, 

without the three-body potential (Z *■ 0) , the diamond cubic and simple 

cubic structures are not stable. However, increasing the three-body 

Intensity parameter, Z , to values greater than Z 0.5, one can 

stabilize the diamond cubic structure. For the simple cubic structure, the 

★ 

elastic constants satisfy all the stability requirements for Z values 

A 

in the range 0.45 < Z < 0.73. 
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Using eqs. (1) and (2) and Figs. 2 , both the thermodynamic stability 

condition (minimum value of and the mechanical stability condition 

* 

may be compared so that the domain of Z providing absolute stability 
for each cubic structure may be determined. These results are shown In 
Flg< 3. From Fig. 3, we note that the FCC structure Is absolutely stable 

ik 

In the range 0 < Z <0.45, the BCC structure has no region of absolute 

stability, the simple cubic structure Is absolutely stable In the range 

0.45 < z'*' < 0.57 and the diamond cubic structure Is absolutely stable In 

the range Z > 0.57. Of course, these ranges of absolute stability 

depend on the choice of crystal structures used In the comparison and upon 

the type of two**body and three**body potentials used in the potential energy 

function. The use of different exponents (m and n) In the two-body 

*(7) 

potential are known to Influence the niagnlcude of but their 

variation has not been analyzed In the present investigation. 

CONCLUSIONS 

Atomistic models based on two-body-only Interactions (Z « 0) cannot 
be used to simulate the various properties of open structures like the 
diamond cubic structure because this structure Is not mechanically stable 
under such conditions. Multi-body Interactions, or at least three-body 
interactions, are absolutely necessary to provide both thermodynamic and 
mechanical stability for such open structures. The domain of absolute 
stability for each crystal structure can be found for a chosen potential 
energy function and chosen list of comparison structures. 
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Table 1 


* 

Lattice Sums of the Energy for Mle and Axilrod-Tellor 
potentials (m ■ 12 and n ■ 6) 



FCC 

BCG 

* 

DIA 

SC 

B 

14.44B1 

12.2495 

5.1153 

8.3994 

A 

12.1319 

9.1141 

4.0389 

1 

6.2021 

1 

19.1697 

14.7719 

1.6647 

6.6138 


Diamond cubic structure. 
Simple cubic structure. 
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Table 2 


Lattice Sums of Elastic Constants for FCC, BCC, Simple Cubic (SC) 
and Diamond Cubic (DIA) Structures (o - 12, n - 6). 


Structure 

^12^^1122^ 

^44^^2323^ 

A (others) 

pq 

®11^®1111^ 

®12^®1122^ 

®44^®2323^ 

B_^(others) 

pq 

'^11^’^lllP 

*^12^^1122^ 

^44^^2323^ 
T^_( others) 

pq 


8.1510 

10.1074 

274.7199 


4.0124 

4.4705 

172.4297 

FCC 


4.0124 

4.4705 

57.0059 


0.0 

0.0 

0.0 


5.0125 

8.1746 

222.0875 


3.5698 

4.0279 

129.8891 

BCC 


3.5698 

4.0279 

40.4850 


0.0 

0.0 

0.0 


8.1337 

9.6340 

96.4657 


0.0679 

0.7563 

59.9206 

SC 


0.0679 

0.7563 

20.9476 


0.0 

0.0 

0.0 


1.8050 

2.6088 

26.0069 


1.7899 

2.0955 

14.1703 

DIA 


1.7899 

2.0955 

9.2187 


0.0 

0.0 

0.0 
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APPENDIX I 






Lattice Sua Fonailac for Elastic Constanta Dslng Hie Potential 
and Azllrod-Tcller Potential (three-bodj) Teras 

Defining our potential, ee the sum of a two~body, $ 2 * 
body, tern, the two-body part can be written as 


N 


$ « ii£ V /JL- f_£) 
2 2 jt, 'm-n ^2^ 


m/2 


iwl 


r^ n/2 

^ f-2.1 I 

m-n ^2^ ' * 

*^1 


Let ue define 


x„.m/2 


«n n/2 


*0-r^ -5^(5^) . 


then we have 


Ne S 


Ne mn r >m 
2 2(m-n) ‘^2 


r m 


r n 


{(, + 1)(^ ) Aj(py5 “ (9 *** DCrf ^ 


oPy6' 


where the two-body lattice sums ®nd are given by 


N 


apyS 


j m 4x? xf xT X? 

r (/-) ^ 

1-2 -'ll r,. 
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(two-body) 


a three- 


(1-1) 


( 1 - 2 ) 


(I-3a) 


(I-3b) 


(I-4a) 


! 




’apv6 


. n 4x® xT X? 

I - - 4“ — 


N 


i-2 ‘11 


11 


(I-4b) 


and 2xJxJ - BX^/&T)p^. 

The three-body part of the potential, given by 


- 


irhll 

^ 1 j 
i*tj 


N N 1 + 3 COB e. COB 0. COB 0, 

f 1 , J L 


<'11 'U '«’■ 


(I-5a) 


N N 


■fr^ini- 


1 


,2.2 2 w 2 . 2 2 ,, 2 . 2 2 . 

+ 2 ^11*^^ 1.1 ~*^H ^ ^ 11~^*^1.1 1.1 ^ ^ 1.1 li “ 1 1 ^ 

j . V J 8 , 2 . 2 2 v5 ' 

(I-5b) 


from Fig. 1. 


2 2 

Let us define X " r^^, Y ■ r^ and Z 


r^j and as 


. __1 . 3 (X+Y-2) (X-f-2-Y) (Y+Z-X) 

*3 ■ <xyz)3/2 8 (XYZ)8/<< 


<I-6) 


Then, the three-body lattice sums in the elastic constants can be written 
as 


^o 0 y 8 



N N 

I I (e^ ^3 e2> 

1 J ^ ^ ^ 

i\l 


(I- 7a) 


where e^ and 62 are the column vectors, Is the (3 x 3) matrix 
and these have the following elements 
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•T / SX 

" te’ 


^3- 



~T 

f ax BY 

a: 

, 

' 

®2 " 


’ an 

a\ 

a^3 

a^B 3 " 


ax"* 

BXBY 

am 


9%3 

a^4>3 

a%3 


BYBX 

a^ 

BYBZ 

• 

a 24,3 

a '‘ 4)3 

th 


azax 

aZBY 

az^ 



<I-7b) 


The elements In eq. (I-7b) can also be written using atomic coordinates: 


— 1. 2x^ ■ 2x** x^ * «■ 2(x ®*-x^(x^“X^) (I-8a) 


'op 


dp 


and 


ax^ 


2 2 

a 4>, 3 4*, 

- £<X,Y,Z) ; - f(Y,Z,X) ; 

air az 


f(z,x,Y) , 


2 2 
a^*3 B^(t>3 


2 2 

a *3 a^(i .3 


axaY " BYax " s^^c.y.z) ; " 5zax “ 8 (z»x,y) , 


2 2 
a^<t>3 a^d.3 

BYaz 


azBY " » 


(l-8b) 


where £(X,Y,Z) and g(X,Y,Z) are de£lned as 
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f 


f(X,Y,Z) 


^^T79 -4 + ^77 [2(V+Z-X> - 2CX+Z-Y) - 2CX+Y-Z) 

4(XYZ)‘^'^ X 8(XYZ)^'^ 


+«! {(X+Y-Z)(X+Z~Y) - (X+Z-Y)(Y+Z-X) - (X+Y-Z) (Y+Z-X) } 


+ “^ {(X+Y-Z)(X+Z-Y)(Y+Z-X)}] , (I- 8 c) 

4X^ 


and 

g(X,Y,Z) 

+ I (X+Y-Z)(Y+Z-X)(| - 7 ) + I (X+Y-ZKX+Z-YK- 7 + 7 ) 

+ 11^ (X+Y“Z)(X+Z-Y)(Y+Z-X) ] . (l- 8 d) 


4 + 


4(XYZ)^'^^ ^ 8 (XYZ) 


[2(X+Y-Z) + I (X+Z~Y)(Y+Z-X)(- 7 - 7 ) 


Using eqa. (7) and ( 8 ), the three-body lattice sums for the various crystal 
structures were calculated and listed in Table 2. 
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FICUBE CAPTIONS 

* 

1. Atom geometry needed In the evaluation of strain parameters and the 
lattice sums of the elastic constants* 

2. Variation of the elastic conatanta Cj 2 » “ 1^121 

(Cjj + 2 Cj 2)/3 *1® ® function of Z* for (a) FCC, <b) BCC, (c) 

Diamond Cubic, and (d) Simple Cubic structures* 

ik 

3. Reduced energy » ^ » for the FCC| BCC, Diamond Cubic (DIA) and Simple 

A 

Cubic (SC) structures as a function of Z , plus mechanical stability 
regions for FCC, fiCC, DIA and SC structures based on their elastic 
constants* 
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FIGURE 2c 
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\/ / /\ and r I indicate the stable and 
unstable regions against elastic deformation, 

respectively 


FIGURE 3 


AA 






Chapter III 


CALCULATION OF POTENTIAL ENERGY PARAMETRES 
FOR THE SILICON - CARBON SYSTEMS 



Introduction 


In simulation calculations, the functional simplicity as well as the reliability 
of a potential energy expression are both important considerations. In general, 
in such modeling investigations based on simple two-body interactions, the close- 
packed structures (hep or fee) are the only energetically stable configurations. Less 
dense structures (such as diamond cubic), on the other hand, are not stable with a 
potential energy function expressed as the sum of pair interactions only. This fact 
has been demonstrated in a recent parametric study on the stability of the diamond 
cubic crystal [l]. It has been determined that multi-body interactions are absolutely 
necessary; therefore, a potential energy function comprising two- and three-body 
terms (at least) must be employed for an adequate description of the diamond cubic 
structure. 

In the present study, specific parameters of the potential energy functions for 
a variety of Stp, C, and Sij^Cy species are calculated and the applicability of the 
potential function, for reproducing various energetic and structural quantities was 
tested. 

Analytical Procedure 


.K. Elemental Systems 

For a system of N particles, in general, the total potential energy may be 
expanded as |2): 

^ ” ^2 "*■ ^3 “^ * * * 


4y 

' X j k 

where, u(f,,fj) and represent the two- and three-body interactions, 

respectively, The position of the i’th particle is denoted by r,-. 

In this study, the two-body part is represented by a Mie-type potential: 
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u{r,;) = 


(m-n) Pr 




tj 


i; 


( 2 ) 


where, r,/ = |ri - fjj; r„ represents the equilibrium distance and e denotes the 
two-body energy at r,’y = to. The exponents m and n account for the repulsive and 
attractive terms, respectively, For the three-body term, the Axilrod-Teller equation 
is taken into consideration: 




Z(1 -f 3 cos Oj cos Oj cos 

(r,; • 


( 3 ) 


where, 6i,0j,6k and r,-j , r,*fc, ryjk represent the angles and the sides of the triangle 
formed by the three particles i, j and k, respectively. The three-body intensity pa- 
rameter is denoted by Z. Interactions coming from higher body terms are neglected 
in this calculation. Considering the lattice sum formalism, next, the total potential 
energy function for a crystal is obtained by combining equations 1 through 3: 


4» = 


Nc 


nAC-fr-mBC-fr 


■^NZTkCj? 


(4a) 


2{m - n.) 

where, d represents the nearest neighbor distance and the lattice sums are given by: 



(46) 

(4c) 


and 



(1 -r 3 cos di cos 9j cos Ok)d^ 
(r«j • Uk • ryfc)3 


(4d) 


Numerical values for the lattice sums A, B and T/t for the diamond cubic crystal 
are given in reference [1). 


B. Binary Species 

For a binary system, equation la is replaced by 

^ 53 S *^3(0:4 (3) 

Q=l/?=1 astl/3=n=l 
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W W4. ,4. MKt- : » . 


where, 0a and 03 are still given by eqs. lb and Ic but now deal with the interactions 
between the species a and 0 and '){') =a or 0), For a binary crystal, we have 


02(cr,/?) = 


(6a) 

where 

N(i , 

(66) 



(6c) 

and 

fl^3(a,/?,Tf) = 

(7a) 

with 

1 ^ ^ (1 + 3 cos Of cos 0j cos 02)d^ 

- 3! ^ ^ ‘ r,k(a,y) • 

’ }ik 

(76) 


Subscripts t, j and k are the running indices, with similar definitions as the variables 
in eq. 4. 


Pertinent lattice sums for several single component and binary systems are 
given in Tables 1 and 2 respectively. The following relationships hold for the two- 
component lattice sums 

An = Az2 i An = A2\ ; .Bn = B22 ; B12 = B21 ; 

Tin = T* 3 , T112 = T121 = T212 - T221 ; T122 = T211 ( 8 ). 


C. .f'arameter Evaluation 

In order for the potential energy function given by equations 1 through 3 to 
be used in a simulation calculation for a specific system, the parameters (c, r„ and 
Z) must be first evaluated. In accordance with our earlier study, on the stability of 
the diamond cubic structure, the values of the exponents m and n were taken as 12 
and 6, respectively. The evaluation process is basically a simple fitting procedure, 
however, due to the nonlinear nature of the potential function It often becomes 
quite cumbersome. Throughout this study only experimental data were employed 
in the evaluation of the parameters. Furthermore, one should remember that the 
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parameters, by definition, are independent of N and of the geometrical state of the 
system (i.c., independent of the particles’ positions), They depend solely on the 
atomic species involved in the interaction. 

In this study, parameters for pure Si and C as well as for the SIC system were 
evaluated, In each case, the basic experimental data for the crystalline bulk and 
for smalt clusters were taken into consideration. The data which are used in the 
calculation of the parameters are tabulated in Table 3, The evaluation procedure for 
the pure systems will be first outlined. The parameters for the homon'Jclear cases 
(Si and C) were calculated considering three basic restrictions, (i) The equilibrium 
criterion; at the static iimit we have d^/drq = 0; where, in the case of a crystal, 
rq is the nearest neighbor distance and is directly related to the density. For small 
clusters, however, r, represent internuclear distances between the parti-^les. (ii) The 
stability region for the diamond cubic structure: according to results reported in 
reference [l] the parameter Z should satisfy: 0,55 < Z” < 0.8, where, Z' is the 
reduced three-body intensity parameter defined as: 2^' = Z/cro®. (iii)Stable trimer 
configuration: the energetically most favorable trimer configuration for a potential 
energy function defined by equations 1-3 has been investigated In reference [3]. 
Accordingly, if the trimer is linear, 0.68 < Z’ and, if it is (equilateral) triangular, 
Z' < 0.68 must be satisfied. 

The parameters for the homonuclear cases (Si and C) were selected in such 
a way that calculated energies for the small clusters (5i3, C3, 5t2 and C2) and 
the crystalline states (silicon and diamond) exhibit a best fit to the corresponding 
experimental data. For this purpose the energies were calculated according to equa- 
tions 1-3 using a series of parameter sets which obey the restrictions cited above. 
Numerical values of the best fitted parameters are given in Table 4. 

In the case of the parameters for the SiC system, experimental data for 5tC, 
512C, SiCi molecules and for the crystalline ff-SiC were used. The data consisted 
of the cohesive energy and the nearest neighbor distances for the bulk, and the 
bond energies along with the bond lengths for the small clusters. (See Table 3). 
Due to the binary nature of this system, in addition to the homonuclear interaction 
parameters, cross-parameters such as csi,Ct fo{si,c)i 2si,si,c and Zsi,c,c needed 
to be determined. 


In the evaluation of the cross-parameters for the 0-SiC system, we also consid- 
ered the following stability restrictions, (i) According to experimental data, St 2 C 
is a linear symmetrical molecule. Our present potential energy function satisfies 
this condition only if > 0,41 where, 


2 '- . . _ Z^si,si,C) 

(St,S:,C) (€(s,'^5,)C(5,-^C-)£(5,-^<-j)>/3(ro(5,-^5,)ro(5,-,c)ro(S,-,C})® 


(9) 


(ii) On the other hand, experiments indicate that SiC 2 is also linear but asymmet- 


49 



I 


rical and that it can only be stable If Z’[S\,c^c) > 0*58 with, 

^{Si,C,C) 

(£(Si,C)«(5i,C)t(C’,C))^'^^(»‘o(5i,r)^o(S«’,r)ro(C',C’))® 

(ill) For the crystalline 0-SiC as well as for the 5tC, ifftjC and iffiCa molecules, 
the stability criterion which is given by = 0 should also be satisfied. 

We employed the parameters for SI and C obtained above, and calculated 
energies for the crystalline and a-5*C, .^i 2 C and StC^ molecules for a 

series of cross-parameters acts. The best fitted cross-parameters with the above 
restrictions were chosen and tabulated in Table 4 along with the parameters for Si 
and C. 



Discussion 

In order to analyze how well the parameters can now reproduce experimental 
data, we plotted the experimental energies, versus the calculated values in figure 
1. It Is important to remember that, in this graph for each case, we used the same 
potential function associated with the same set of parameters given in table 4, The 
dotted line indicates the ideal positions (i.e., the 45‘’ line). For such a large variety 
of systems (ranging from crystalline phases all the way down to isolated molecular 
species) the agreement obtained between the experimental and calculated energy 
values can be considered quite good. The plot also contains the energies for C 4 , 
Cs and crystalline graphite which were not included in the parameter evaluation 
procedure. While the energy value for graphite displays some (~ 30%) deviation, 
calculated energies for C 4 and Cs are in very good agreement with the experimental 
values. Another important point related to figure 1 is the slight segregation of the 
energy values for different species. All the calculated energies for the bulk cystalline 
cases (namely, for silicon, diamond, graphite, a-silicon carbide(4H-III) and /^-silicon 
carbide) were found to be somewhat lower than the corresponding experimental 
values, while all the small molecular species produced slightly higher energies. In 
general, the tendency for this sort of segregation in the case of a potential energy 
function based on two-body interactions only, is much more pronounced [4] than in 
the present case. This situation may be closely related to the effect of three-body 
Interactions which are properly accounted for in our calculations via equation 3. 
The energies for o-SiC and /?-SiC were found to be very close to each other (see 
Table 3), and they were considerably lower than the energies calculated for other 
structures (Table 2). 
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Figure Caption 


Figure 1 Comparison of the calculated cohesive energies with the experimental 
values. 
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Tabic 1 


Lattice Sums for the Mie and Axllrod-Teller Potentials for Various Structures. 



FCC 

BCC 

DIA 

GRAPH 

/?-tin 

B 

14.4481 

12.2495 

5.1153 

3.3895 

8.2864 

A 

12.1319 

9.1141 

4.0389 

3.0092 

5.4654 ' 

Tk 

19.1697 

14.7719 

1.6647 

Q.lOlO 

7.0706 


Table 2 

Lattice Sums for two-component systems. 


NaCl 

CsCI 

j ZnS 

BN 


structure 

structure 

• structure 

structure 


”0.1896 

1.1038 

• 0.0337 

0.0084 

i 0.0342 

6.0126 

8.0103 

. 4.0045 

3.0009 

: 4.0595 

1.8003 

3.5357 

, 0.7596 

0.2821 

; 0.7657 

6.5888 

8.6996 

! 4.3516 

3.1163 

' 4.3893 

2,5272 

5.4138 

0.6928 

0.1057 

' 0.7013 

11.5117 

25.8716 

; 2.8595 

0.1574 

i 2.8967 j 

5.7544 

12.9374 

1.4298 

0.0787 

‘ 1.4558 1 
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Table 3 

Experimental and calculated energies and bond distances. 


Species 


S>3(p) 

linear 


Si(,) , 

diamond cubic 


Cj 


Itot’r 


,?4(p) 

linear 


C6(ff) 

linear 


diamond cubic 


graphite 


SiC(„) 


symmetric linear 


CsSijp) ^ 

unsymmetne linear 


/?-SiC(,) 

zinc-blende 


Q-SiC(a) 

4H-1II 


Bond distances or nearest 
neighbor distance (A 


Experiment Calculation 


2.246< 2,295 


2.25' 2.267 


2.3516 


1.4806 


1.458 


n = 1,44 
f2 ~~ 1.455 


n = 1,435 
T2 = 1.455 


1.5445 


1.4619 



1.4180'® 



1.885'®^ 


C-C- 1.463 
C-Si=1.7ll 


1.8115 


1.8135 


Cohesive energy 
(cV/molecule) 


Experiment 


-3.208 ±0.216 I -2.817 


-7.407 ± 0.563 


•4.662 ± 0.043 


- 6.2 


-14.033 


-19.686 


-27.014 


saiKctfVi 


- 0,007 






- 0.309 


-7.394 


-7.415 


- 0.022 


- 0.022 


-13.128 


- 0.300 


■12.835 


• 12 * 817:^:,30 



-5.568 


-5.437 


-12.362 


-26,70 


-9,402 


-9.807 


-3.894 




















































Tabic 4 

Potential Energy Parameters. 


Two-body 

Parameters (Si-Si) (C-C) (Si-C) 

c (eV) 2.817 5.437 3.895 

ro (A) 2.2951 1.4806 1.74 


Three-body 

intensity 

parameter (Si-Si-Si) (C-C-C) (Si-Si-C) (Si-C-C) 


Z (eV-A®) 3484.0 167.3 698.2 261.8 
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